Topological Order from Geometric Entanglement 
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Here we investigate the connection between topological order and the geometric entanglement, 
as measured by the logarithm of the overlap between a given state and its closest product state 
of blocks. We do this for a variety of topologically-ordered systems such as the toric code, double 
semion, color code, and quantum double models. As happens for the entanglement entropy, we 
find that for large block sizes the geometric entanglement is, up to possible subleading corrections, 
the sum of two contributions: a non-universal bulk contribution obeying a boundary law times the 
number of blocks, and a universal contribution quantifying the underlying pattern of long-range 
entanglement of the topologically-ordered state. This topological contribution is also present in the 
case of single-spin blocks, and constitutes a novel characterization of topological order based on a 
multipartite entanglement measure. In particular, we see that the topological term for the 2D color 
code is twice as much as the one for the toric code, in accordance with recent renormalization group 
arguments [H. Bombin, G. Duclos-Cianci, D. Poulin, New J. Phys. 14 (2012) 073048]. Motivated 
by these results, we also derive a general formalism to obtain upper- and lower-bounds to the 
geometric entanglement of states with a non-Abelian group symmetry, and which we explicitely use 
to analyze quantum double models. Furthermore, we also provide an analysis of the robustness of 
the topological contribution in terms of renormalization and perturbation theory arguments. Some 
of the results in this paper rely on the ability to disentangle single sites from the quantum state, 
which is always possible for the systems that we consider. Finally, in the Appendix we relate our 
results to the behavior of the relative entropy of entanglement in topologically-ordered systems. 



I. INTRODUCTION 

Topological order (TO) [ ' ] is an example of new physics 
beyond Landau's symmetry-breaking paradigm of phase 
transitions. Systems exhibiting this new kind of order are 
linked to concepts of the deepest physical interest, e.g. 
quasiparticle anyonic statistics and topological quantum 
computation [ ]. Importantly, TO finds a realization in 
terms of topological quantum field Theories ['], which 
are the low-energy limit of quantum lattice models such 
as the toric code and quantum double models [4] , as well 
as string-net models ["]. 

A remarkable property about TO is that it influences 
the long-range entanglement in the wave function of the 
system. For instance, as proven for systems in two spatial 
dimensions (2D) [6, 7] , the entanglement entropy S" of a 
region, such as A in Fig. l.a, of boundary size L ^ 1 
obeys the law 



S^So-S^ + OiL-"), 



(1) 



where 5*0 oc L is a non-universal term (the so-called 
"boundary law", or "area law"), v some exponent, and S-y 
is a universal long-distance contribution: the topological 
entanglement entropy. S-y is non-zero for systems with 
TO, e.g. Sy — I for systems in the topological phase of 
the toric code. In general one has that 



(2) 




where {da} are the so-called quantum dimensions of the 
associated anyon model. More recently, similar universal 
contributions have also been found for other bipartite 
entanglement measures such as the mutual information 
and the Renyi entropy [8, 9]. 

The main purpose of this paper is to investigate, in 
the context of systems with TO, a global measure of 
entanglement which captures truly non-bipartite corre- 
lations in the system: the geometric entanglement (GE) 
[10, 11], which we shall denote by Eq. This measure in- 
tuitively characterizes how well an entangled state can be 
approximated by a mean-field state, i.e., a product state. 
Can such a mean-field, global measure of entanglement 
be able to reveal any property about topological order? 
As we shall see in this paper, this is indeed the case. 

Most of our study is focused on exactly-solvable mod- 
els corresponding to fixed-points of the renormalization 
group (RG). We deal most prominently with the toric 
code model [4], but shall also discuss the double semion 
model [•")], topological color codes [12], as well as quantum 
double models [4] . Each one of these models can be recast 
as some particular instance of a string-net (Levin- Wen) 
model. Because of this, they correspond (by definition) 
to RG fixed-points and, therefore, they are representative 
of their respective topological phases. Our main result 
is that for these various fixed point models, the GE of 
ground states obeys 

Eg = Eq - E^, (3) 

with E^ being identical to S^ and Eq some non-universal 
term. The exact form of Eq depends on short-distance 



details such as the size of the blocks (i.e. the number of 
spins) considered for the closest product state optimiza- 
tion. As we shall see, for blocks of boundary size L, Eq 
will be proportional to L and also to the number of blocks 
n\, (more specifically, Eq oc ni,L). We call this behavior, 
which depends on short-distance correlations, a "bound- 
ary law" for the GE, since it corresponds to a boundary 
law term for each one of the blocks. Moreover, E^ being 
identical to Sj is an indication of the topological origin 
of the term E^ . 

Additionally, we argue that the above law is also valid 
up to subleading corrections away from the RG fixed 
points, e.g. for the toric code model under the influ- 
ence of non-relevant (and weak) perturbations such as 
magnetic fields. In this case the geometric entanglement 
of blocks with boundary size L ^ 1 seems to obey the 
law 



bipartite 



Eg = E^-E^ + 0{L-''), 



(4) 



where again Eq cx n^L is a non-universal term, v' some 
exponent, and E^ the topological contribution. From 
here on, we shall refer to the term E^ as the topological 
geometric entanglement. 

It is worth mentioning now a number of remarks about 
our results. First, they provide the first evidence that 
topological order can actually be read and assessed from 
the geometric entanglement, which is a multipartite (in- 
stead of bipartite) measure of entanglement. Thus, this 
may be considered as a new way of characterizing topo- 
logical order in quantum states. Second, our results arc 
also the first explicit example of a boundary law behavior 
for the GE in ground states of 2D quantum many-body 
systems. Third, we will see that the topological term for 
topological color codes is twice the one for the toric code. 
This is consistent with the recent result that the former 
model is equivalent to two copies of the latter under some 
RG transformation [L'i]. 

Importantly, there may be potential advantages in con- 
sidering the GE to characterize topological order instead 
of other quantities. For instance, one of them is that 
given a (possibly approximate) description of the ground 
state by a Tensor Network such as a Projected Entan- 
gled Pair State (PEPS) [14] or a Multiscale Entanglement 
Renormalization Ansatz (MERA) [15], its computation is 
not too costly. This is because, in the context of Tensor 
Network methods, the calculation mainly involves over- 
laps with product states and optimizations over these, 
which can be computed quite efficiently. But moreover, 
the GE can also be probed experimentally with current 
technology using e.g. Nuclear Magnetic Resonance via 
single spin measurement on all nuclear spins [Ifi], or ul- 
tracold atoms in optical lattices via similar techniques to 
those explained in Ref. [j V]. 

The structure of the paper is as follows. First, we 
briefly review the basics of the geometric entanglement 
in Sec. II. After this, we deal in Sec. Ill with the toric 
code model. Since this is the simplest model displaying 
non-trivial topological order, we spend quite some time 
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FIG. 1: (color online) Two different types of partitions, (a) 
Bipartition with region A and region B. Region A has linear 
size L. This is usually the partition used in evaluations of the 
entanglement entropy and other bipartite entanglement mea- 
sures, (b) Multipartition with many blocks, each one with 
linear size L. This is the partition used in the calculations of 
the geometric entanglement. If each block contains only one 
spin (akin to individual degree of freedom) , then the product 
state is a completely separable state and the geometric entan- 
glement is defined with respect to individual spins. However, 
each block can also contain multiple spins, and in this case the 
geometric entanglement is defined with respect to the closest 
product state of the individual blocks. In this latter scenario, 
we can inquire how the entanglement depends on the block 
boundary size L and, in particular, whether there is some- 
thing like a boundary law (i.e. linear scaling with L) and 
whether there is a topological correction. 



explaining many of its properties, as well as derivations 
of the geometric entanglement of spins and blocks both 
for the square and honeycomb lattices. In Sec. IV we 
explain a number of important points that need to be 
considered in order to generalize the results obtained for 
the toric code to other models. Then, the double semion 
model is analyzed in Sec. V, and topological color codes 
are addressed in Sec. VI. For color codes we derive results 
using two independent methods: a bound approach (sim- 
ilar to the one used for the other models), and a direct 
approach valid for Calderbank-Shor-Steane (CSS) self- 
orthogonal codes. Depending on the particular setting we 
will see that one approach may provide some advantages 
over the other. After this, we develop in Sec. VII a gen- 
eral formalism for the bound approach, valid for states 
with non-Abelian symmetries, and apply this formalism 
to quantum double models in Sec. VIII. The situation be- 
yond RG flxed points is considered in Sec. IX, where the 
robustness under perturbations of the topological contri- 
bution is assessed by RG and perturbation theory argu- 
ments. Finally, we present our conclusions and possible 
future directions in Sec. X, and also describe in the Ap- 
pendix how our results are related to those of the relative 
entropy of entanglement for topologically-ordered states. 



II. GEOMETRIC ENTANGLEMENT 

To begin with, let us recall some basic notions about 
the geometric entanglement. Consider an m-partite nor- 



malized pure state |*) e H = 0^1^^, where "HW is 
the Hilbert space of party i. For instance, in a system of 
n spins each party could be a single spin, so that m — n. 
But each party could also be a set of spins, either con- 
tiguous (a block [is]) or not. 

We wish now to determine how well state |^) can be 
approximated by an unentangled (normalized) state of 
the parties, |<I>) = (8)™i |(/)t'l). The proximity of |^) to 
|<I>) is captured by their overlap. The entanglement of 
1^) is thus revealed by the maximal overlap [10], 



Amax(*) 



max|(<f>|*) 



(5) 



The larger A^ax is, the less entangled is |^). Therefore, 
we quantify the entanglement of |^) via the quantity 



i?GW = -l0g2ALxW: 



(6) 



where we have taken the base-2 logarithm, and which 
gives zero for unentangled states. £'g(^) is called geo- 
metric entanglement (GE). This quantity has been stud- 
ied in a variety of contexts, including critical systems 
and quantum phase transitions [18, 19], quantification 
of entanglement as a resource for quantum computation 
[20], local state discrimination ["Jl], and has been recently 
measured in NMR experiments [Ki]. 

To make a connection between the GE and other mea- 
sures of entanglement, let us consider the case of just 
two sets of spins. In this case £'g(^) coincides with the 
so-called single-copy entanglement between the two sets. 



Ei{-^)^-\og^i'i{p), 



(7) 



with vi (p) the largest eigenvalue of the reduced density 
matrix p of either set [22]. As is well known, this also 
coincides with the a-Renyi entropy 



Sa 



1 



(l-«) 



log(tr(p")). 



(8) 



for a — )■ oo. 

Unlike other measures of entanglement , the GE offers 
a lot of flexibility to study multipartite quantum corre- 
lations in spin systems. For instance, one can choose 
each party to be a single spin, but one can also choose 
blocks of increasing boundary length L [18] (see Fig. l.b 
for an example). In fact, studying how the GE changes 
with L provides valuable information about how close 
the system is to a product state under coarse-graining 
RG transformations. 



III. TORIC CODE MODEL 

As described in the introduction, here we considered a 
variety of models corresponding to different topological 
phases. The simplest of these models is the toric code [4]. 
This model is, in fact, the simplest example of a topo- 
logically non-trivial 2D system, and our results for this 



model provide also the grounds for the rest of the mod- 
els that we shall consider (namely double semion, color 
codes, and quantum double models). Given its relevance 
for the rest of the paper, we provide now a short intro- 
duction to some of the key aspects of this model. As 
we shall see, this will be very useful for the forthcoming 
calculations. 

Mathematically speaking, the toric code is the RG 
fixed point of the topological phase of a Z2 gauge the- 
ory. The model is equivalent under local transformations 
(or local moves, or local disentanglers) to the Levin- Wen 
string model on a honeycomb lattice [5], which is by def- 
inition a RG fixed point. A more precise derivation of 
this property will be provided later. 

To define the toric code, we consider a lattice S on a 
torus. In this paper we will deal with the toric code in 
the square and honeycomb lattices, yet the model can be 
defined on arbitrary lattices. Other Riemann surfaces of 
genus could also be considered easily without changing 
our conclusions. There are spin-1/2 (qubits) degrees of 
freedom attached to each link in lattice S. The model 
is described in terms of stars and plaquettes. A star "s" 
is a set of links sharing a common vertex. A plaquette 
"p" is an elementary face on the lattice E. For any star 
s and plaquette p, we consider the star operators Ag and 
plaquette operators B^ defined as 



A_, 






■m 



3&P 



(9) 



Jj] 



where aa is the a-th Pauli matrix at link j of the lattice. 
Let us call respectively Ug, Up and n the number of stars 
(vertices), plaquettes (faces) and links (sites) in lattice 
S. Importantly, star and plaquette operators satisfy the 
global constraint 



n^-n^p 



I. 



(10) 



Therefore, there are ris — 1 independent star operators 
and Up— 1 independent plaquette operators. 



With the definitions above, 
model reads 



the Hamiltonian of the 



Htc 



Y^Ag-Y^Bp. 



(11) 



This Hamiltonian is frustration free (i.e. all the terms 
in the above sum conmute with each other), and can 
be diagonalized exactly as explained in Ref. [ ]. The 
ground level is 4-fold degenerate on a torus (40-fold for 
a Riemann surface of genus g). This degeneracy depends 
on the underlying topology of lattice S, and is already 
by itself a signature of topological order. Moreover, the 
ground level is a stabilized space of Gs , the group of all 
the possible products of independent star operators, of 
size I^J = 2("=-i). 



A. Ground states 

In order to build a basis for the ground level subspace, 
let us consider a closed curve 7 running on the links of the 
dual lattice S*. We define its associated loop operator 

^a; [7] = niG7 "x , where j G 7 are the links in E crossed 
by the curve 7 connecting the centers of the plaquettes. 
It is not difficult to see that the group Gs can also be 
understood as the group generated by all the possible 
contractible loop operators [ ]. Let also 71 and 72 be 
the two non-contractible loops on a torus, and define the 
associated string operators wi^2 = I^2;[7i,2]- We call |0) 
and |1) the eigenstates of cTz respectively with +1 and 
— 1 eigenvalue. With this assumptions, the ground level 
subspace then reads C — span{\i, j), i,j — 0, 1}, where 




\hJ) 



V\Gs 



see. 



Wr 



|0) 



0n 



(12) 



It is easy to check that the four vectors \i,j) are orthonor- 
mal and stabilized by Qs, that is, g\i,j) — \i,j) Vg € Gs 
and yi,j. These four states form a possible basis of the 
ground level subspace of the toric code model on a torus. 



B. Excited states 

Excited states of the toric code can be constructed by 
locally applying Pauli operators a^^o-y, a^ on the ground 
states \i,j). Pauli operators az create pairs of decon- 
fined charge-anticharge quasiparticle excitations, ax cre- 
ate pairs of deconfined flux-antiflux quasiparticles, and 
ay creates a flux-antiflux and a charge-anticharge pairs. 
These operators can be applied to several sites of the lat- 
tice, thus creating a quasiparticle pattern that defines the 
excited state. The excited states of the model are then 
labeled by a set of quantum numbers, \(j),c,i,j), where 
(j) and c are patterns denoting the position of flux-type 
and charge- type excitations respectively, and «,j label 
the ground state \i,j) that was excited. In this paper 
we will focus on the entanglement properties of states 
l<A,c,i,j). 



C. Disentangling the toric code 

A key property of the toric code model, which is fun- 
damental for some of the derivations in this paper, is that 
its ground state |0,0) = |z = 0, j = 0) can be created by 
a quantum circuit that applies a sequence of ControUed- 
NOT (CNOT) unitary operations over an initial sepa- 
rable state of all the qubits [2-'>]. This means that it is 
actually possible to disentangle qubits from the ground 
state of the system simply by reversing the action of these 
CNOTs. Moreover, these CNOT operations are local, i.e. 
they do not span over delocalized sites in the lattice, but 
rather act over pairs of not-too-far-neighbouring qubits 
(e.g. nearest- and next-to-nearest-neighbours) . This was 



FIG. 2: (color online) (a,b) CNOT operations that disentan- 
gle qubits from the system, thus removing their links from 
the lattice. In the diagram, the arrows go from controlling 
to target qubits; (c) Example of a partition into blocks of a 
given boundary L, e.g. L = 12 here (but all our derivations 
work for arbitrary L). Qubits crossed by the boundary of 
a block are regarded as being inside the block; (d,e) Doing 
CNOT operations locally inside of each block we can remove 
all the stars inside of all the blocks in two steps: first, we ap- 
ply CNOTs as in (a) to get (d) and, second, we apply CNOTs 
as in (b) to get (e). Notice that in (e) each star is made of 
either 4 or 8 qubits. The procedure works similarly for other 
blockings. 



the key observation that allowed to build an exact MERA 
representation of the ground states \i,j) of the model 
[ - ] . The two fundamental disentangling moves are rep- 
resented in Fig 2.a-b, and leave the overall quantum state 
as a product state of the disentangled qubits with the rest 
of the system. What is more, the rest of the system is left 
in the ground state of a toric code model on a deformed 
lattice E, where S is obtained from S by removing the 
links that correspond to the disentangled qubits. This 
property turns out to be of great importance for some of 
our derivations. 



D. GE of the toric code 

We are now in position to study the GE of the toric 
code model. Given that this is a paradigmatic model 
of topological order, we will perform a detailed analysis. 
First, we will provide a number of lemmas and theorems 
that will bring us towards an expression for the GE of 
the toric code in the square lattice, both for spins and 
blocks. Second, we will make a similar study to the case 
of the honeycomb lattice. There we will see that different 
choices of ground states may give rise to different expres- 
sions for the GE of spins, yet the long-range properties 
are the same and equal to those for the square lattice 
after appropriate blocking. This is expected, since both 
lattices can be mapped between each other by means of 



(b) 





FIG. 3: (color online) Bipartite lattices of spins and blocks: 

(a) Sets of spins A and B for Theorem 1 and a 4 x 4 torus; 

(b) Sets of blocks A and B for Theorem 2. 



CNOT disentangling operations (as we will show), and 
therefore correspond to the same RG fixed point. 



1. A couple of Lemmas 

Let us start by considering two Lemmas that will be 
quite useful in our derivations: 

Lemma 1. All the eigenstates \(j),c,i,j) of the toric code 
Hamiltonian have the same entanglement properties. 

Proof. Let us choose a reference state from the ground 
level basis, e.g. |0, 0). First of all, notice that the four 
basis states \i,j) in the ground level are related to |0, 0) 
by local unitary operators, since the string operators wi 
and W2 are tensor products of 2 x 2 identity operators 
and Pauli-a: matrices. Since these operators act locally 
on each spin, they do not change the entanglement of the 
quantum state. Moreover, the excited states \4>,c,i,j) 
are created locally by applying tensor products of 2 x 2 
identity and Pauli-x and z operators to \i,i). The states 
resulting from these operations have then the entangle- 
ment properties of state |i, j), which in turn are the same 
as those of the reference state |0,0). Thus, all states 
|(/),c, i,j) have the same entanglement properties. D 

Lemma 2. Consider an arbitrary set of blocks of qubits 
in lattice S, where each block is regarded as an individual 
party. Then, the entanglement properticsjjf the ground 
state |0, 0) are the same as those of state |0, 0), the ground 
state of the toric code model in the deformed lattice S 
obtained after disentangling as many qubits as possible 
using CNOTs locally inside of each block. 

Proof. Given the ground state |0, 0) in E, we start by 
doing CNOT disentangling operations locally inside of 
each block in order to disentangle as many qubits as 
possible. As a result, state |0,0) transforms into state 

|0'0)di«cntanglcd = l^l) <^ ' ' ' <^ K) ® %Q), whcrC |efe) 

is the quantum state for the k-ih disentangled qubit, 
k = l,...,p (p is the number of disentangled qubits), 
and |0,0) is the i = 0,j = ground state of a toric 
code Hamiltonian in the deformed lattice S. Since all 



CNOT operations are done locally inside of each block, 
states |0,0) and |0, O)^;^^,^^^. j^j^j have the same entangle- 
ment content if we regard the blocks as individual parties. 
What is more, the entanglement in the system is entirely 
due to the qubits in |0, 0), which proves the lemma. D 



2. Square lattice 

We are now in position to consider the entanglement 
properties of the four ground level states \i,j) as well as 
the excited states |^, c, i,j). Our main result here is that, 
for these states, the GE of blocks consists of a boundary 
term plus a topological contribution. Keeping the above 
two lemmas in mind, we now present the following the- 
orem about the geometric entanglement of spins in the 
square-lattice toric code: 

Theorem 1. For the toric code Hamiltonian in a square 
lattice E, the four ground states \i,j) for i,j = 0, 1 and 
also the excited states \(J3,c,i,j) have all the same GE of 
spins and is given by 



Eg = Us- 1, 
where Ug is the number of stars in E. 



(13) 



Proof. Using Lemma 1 we can restrict our attention to 
the GE for the state |0, 0). In this setting, we call compu- 
tational basis the basis of the many-body Hilbert space 
constructed from the tensor products of the {|0), |1)} lo- 
cal basis for every spin, which is an example of product 
basis for the spins. 

Our proof follows from upper and lower bounding the 
quantity Eg{Q,Q) in Eq.(G) for the ground state |0,0). 
First, from the expression in Eq.(33) for the ground 
states we immediately have that the absolute value of 
the overlap with any state of the computational basis is 
l^/^l"^/^, and therefore A^ax > |5s|~^^^- From here, we 
get i?G(0,0) < log2 l^sl) which gives an upper bound. 

Next, to derive a lower bound we use the fact that 
if we group the n spins into two sets A and B, then 
Amax ^ AJnax , whcrc [A : B] means that a partition 
of the system with respect to the two sets A and B is 
considered. Thus, we have that ^g(0, 0) > ^g(0, 0)[^^-^l. 

The trick to find a useful lower bound is to find an 
appropriate choice of sets A and B. In our case, we 
consider e.g. the bipartition shown in Fig. 3. a for even 
X even lattices (other cases can be considered similarly) . 
Then, we use a Lemma by Hamma, lonicioiu and Zanardi 
in Ref. [(>], that the reduced density matrix pA of |0,0) 
for subsystem A satisfies p\ = {\Qs{A)\\Qs{B)\/\Qs\)pA, 
where Qs{A/ B) is the subgroup of Qs acting trivially on 
subsystem A/B. Importantly, for A and B chosen as ex- 
plained above, it happens that Q^ {A) and Qs (B) are triv- 
ial groups consisting of only the identity element. With 
this in mind, we see that the reduced density matrix pA 
has eigenvalues either zero or |t/s |~^ , which is \Qs |-fold de- 
generate. This means that (Amax )^ = l^sl^^ and hence 

^G(0,0)>log2|e.|. 



Combining the two bounds, we get Ec{0,0) — 
log2 \Gs\ = "-S — 1, and from here Eq.(13) for the GE 
for spins follows immediately. D 

Importantly, the techniques used in Theorem 1 can 
also be used to deal with blocks of spins whenever the 
blocks form a bipartite lattice. In general, one may first 
disentangle as many qubits as possible inside the blocks. 
Then the remaining qubits are left in an entangled state 
where the GE is equal to the number of remaining in- 
dependent star operators, which amounts to a boundary 
law term for each block (possibly with a subleading cor- 
rection depending on how the blocks are chosen) plus a 
topological term. As an example of this, let us present 
the following theorem: 

Theorem 2. Given a partition of the square lattice E 
into ni, blocks of boundary size L as indicated in Fig. 2. c 
(where L is measured in number of qubits), then the four 
ground states \i,j) for i,j = 0,1 and also the excited 
states \(f>,c,i,j) of the toric code Hamiltonian have all 
the same GE of blocks and is given by 



Ea = 



ribL 



-1. 



(14) 



Proof. First, notice that Lemma 1 and Lemma pimply 
that we can entirely focus on the ground state |0, 0) of a 
toric code model in the deformed lattice S from Fig. 2.e. 
As shown in Fig. 2.d-e, it is always possible to remove all 
the stars inside of each block in S just by doing CNOTs 
locally inside of each block. As a result, the stars in 
the deformed lattice S correspond to those in S that lay 
among the blocks. It is easy to see that, for n;, blocks of 
boundary L, there are ris = ribL/A of such stars. 

The rest of the proof follows as in Theorem 1, by up- 
per and lower bounding Eq{Q,{)) for the partition with 
respect to blocks. For the upper bound, we use the fact 
that aSLT'^^I > AliiPx'l > \Qs\-^'\ where the first in- 
equality again used the property that if larger blocks are 
considered then the maximum overlap is also larger, and 
Qs is the corresponding group of contractible loop oper- 
ators on E. From here Eq{Q,Q) < logj \Gs\ follows. For 
the lower bound, we divide the system into two sets A 
and B of blocks as indicated in Fig.S.b. Following a sim- 
ilar reasoning as in Theorem 1, it is possible to see again 
that no element g Cz Gs will act trivially on A or B except 
for the identity element. From this point, the rest of the 
proof is simply equivalent to the proof for Theorem 1 . D 



3. Interpretation of the previous Theorems 

Let us now discuss Eqs. (13) and (14) in detail. First, 
let us remind that a variety of works have shown the exis- 
tence of a "boundary law" for the entanglement entropy 
of many 2D systems [25], including models with TO [ , ]. 
It is then remarkable that, according to the first term of 
these equations, the entanglement per block obeys also 



a boundary law, i.e: it is proportional to the size L of 
the boundary of the block. To the best of our knowl- 
edge, these results are the first example of a boundary 
law behaviour for a multipartite (rather than bipartite) 
measure of entanglement in 2D. 

However, the second term in Eqs. (13) and (14) is far 
more intriguing and important. Its existence is caused 
by the global constraint from Eq. (10) on star operators 
which, in turn, allow for the topological degeneracy of 
the ground state. Thus, this term is of topological na- 
ture, universal, and quantifies the pattern of long-range 
entanglement present in topologically ordered states. 

In order to clarify further the meaning of the topolog- 
ical term, let us consider the bipartite case of a block of 
spins of boundary size L and its environment (i.e. the 
rest of spins). In the bipartite case, the geometric entan- 
glement Eg coincides with the single-copy entanglement 
El = — log;/i(/3), with i'i{p) the largest eigenvalue of the 
reduced density matrix p of the block. Since this density 
matrix is proportional to a projector (see Ref. [G] and also 
Ref. [9]), we have that Ei — S, with S the entanglement 
entropy of the bipartition. Thus, for this case we have 
the chain of equalities 



^G 



El 



(15) 



We also know that for a system with TO, the entangle- 
ment entropy satisfies S ^ Sq — Sy -\- 0{L~'^), where So 
is some boundary law term and S^ is the topological en- 
tropy. As explained in Ref. [n], for the toric code model 
the global constraints in Eq.(lO) imply that Sj = 1. 
Thus, with the convention from Eq.(4), we have that the 
topological geometric entanglement is given, in this case, 

by 



E^j 



O-^ . 



(16) 



The topological origin of E^ is thus clear. Notoriously, 
this topological contribution is maintained when promot- 
ing the entanglement measure from the bipartite to the 
multipartite scenario. 



4. Honeycomb lattice 

For the honeycomb lattice, most of the derivations are 
equivalent to the square lattice, yet there is a small dif- 
ference: while the square lattice is self-dual, the honey- 
comb lattice is not. This means, among other things, that 
the toric code in the honeycomb lattice is not self-dual 
with respect to a change of star- and plaquette-operators 
(which is actually the case for the square lattice). In 
practice, this means that different choices of ground state 
basis have different values of the GE of spins. However, 
this is only a short-distance property, since after one RG 
step the toric code in the square lattice is recovered, and 
all the properties of the previous section (including the 
long-distance ones such as E^) hold. In order to show this 
more explicitely, we provide calculations for two specific 
ground states. 



The |0,0) ground state.- We first consider the |0, 0) 



ground state obtained as 
|0,0) = 




0n 



(17) 



which is nothing but the i = 0,j = ground state in the 
notation of Eq.(33). This time, Qs is the group gener- 
ated by the product of all possible star operators in the 
honeycomb lattice. As before, \Gs\ = 2"'="^, where n<, is 
the number of starts (vertices) in the honeycomb lattice. 
For this state, we provide the following observation (not 
a Theorem): 

Observation 1. For the toric code Hamiltonian in a 
honeycomb lattice E, the GE of spins for the ground state 
|0,0) is upper- and lower-bounded by 



Up — I < Eg < ?^p + 1 . 



(18) 



Where Up is the number of plaquettes (faces) in E. More- 
over, numerically it looks like the upper bound is satu- 
rated, so that 



En = n„ 



1. 



(19) 



The above is an Observation, and not a Theorem, since 
Eq.(19) has only been checked numerically. Yet, let us 
prove now the upper and lower-bounds in Eq.(36). First, 
the maximal overlap can be easily upper bounded again 
by Amax > \Gs\^^^^ = 2(i~"=)/2. Nevertheless, unlike in 
the case of the square lattice, this time we can actually 
have a better bound, i.e., 



66^ 






FIG. 4: Bipartition of edges on the honeycomb lattice (or 
equivalently the brickwall lattice) E. We have two sets: A 
(red-grey) and B (blue-white). Edges crossed by the block 
boundaries are regarded as inside the block. 
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16 


24 


8 


9 


4x4 


16 


24 
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4x6 


24 


36 


12 
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TABLE I: GE computed numerically for finite honeycomb 
lattices. We use a, N x N brickwall structure to represent the 
honeycomb lattice with periodic boundary conditions, where 



Us = N^,n= (3/2)A^, and rip 
calculations are accurate with a 10 



n — n. 

6 



(with even A'^). Our 
of relative error. 
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(20) 



which is also smaller than \Qs\ -'-/^ ~ 2^^ ns)/2 ^ging the 

2np for the honeycomb lat- 



fact that n 



Sup and n^ 



tice. This gives the desired upper bound on the GE, 
Eg <np + l. 

To obtain the lower bound on GE (or equivalently up- 
per bound on the maximal overlap) , our strategy is again 
similar to the one for the square lattice. This time we 
consider the bipartition illustrated in Fig. 4. There is a 
star operator in each of the smallest units of the sets A 
(red) and B (blue). It is easy to convince oneself that 
the nontrivial subgroup Gs{A) of Gs is generated by the 
ns/4 star operators residing in all units in A, and sim- 
ilarly for Gs{B). The size of these groups is therefore 
\Gs{A)\ ^ \Gs{B)\ = 2"''/4. The reduced density matrix 
of e.g. set A satisfies p\ = {\Gs{A)\\Gs{B)\/\Gs\)pa, and 
hence the maximal overlap A^^^x is upper bounded by 



A^ < 

max — 



1^.(^)11^.(^3) 



1 



(21) 



which leads to 



1 < Eg, as claimed. Combining this 



bound with the one computed previously, we arrive at 
the result that we mentioned before, namely, Up — I < 
Eg <np-\-l. 

In order to see whether any of these bounds is tight, 
we have carried out a numerical computation for small 
system sizes and confirmed that it is the upper bound 
which is saturated, i.e.. Eg = rip -\- 1. This can be seen 
in Table I. 

The |+,+) ground state.- We now consider the |+,+) 
ground state. By definition, this ground state is obtained 
by acting on the reference state |-|-) " with the elements 
of the group Gp of plaquette operators, namely. 



E 



\^n 



(22) 



=^^1 SGSp 



Notice that, unlike in the square lattice (where there is a 
duality between stars and plaquettes), in the honeycomb 
lattice the |+,+) ground state does not need to have, 
necessarily, the same properties as the |0, 0) ground state 
from the previous section (at least not the same short- 
distance properties). For this state we can actually find 
its GE exactly, and provide the following Theorem: 
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FIG. 5: (color online) (a) Honeycomb lattice (or birckwall 
lattice) E in black lines, and its dual triangular lattice E* in 
dotted red lines, (b) Bipartition of edges in A and B sets on 
the triangular lattice (deformed to fit in to a square lattice). 
For simplicity of the figure, only the spins in the A set are 
shown (in red). The rightmost edges are identified with the 
leftmost edges, and the bottom edges are identified with the 
top edges, since periodic boundary conditions are used. 



Theorem 3. For the toric code Hamiltonian in a hon- 
eycomb lattice E, the GE of spins for the ground state 
|+, +) is given by 



Er 



1, 



(23) 



where rip is the number of plaquettes (faces) in S. 



Proof. We can obtain the GE analytically as follows: 
first, we use the duality that the sate |+,+) in the hon- 
eycomb lattice E is equivalent to the |0,0) ground state 
of the toric code on the dual lattice S*, which is the tri- 
angular lattice. After this, we simply follow the same ap- 
proach as before by finding appropriate upper- and lower- 
bounds. The trivial upper bound is given by rip — 1 < Eq, 
with Up the number of plaquettes in the original honey- 
comb lattice. To find the lower bound, we consider the 
bipartition given in Fig. 5, which is similar to the one that 
we used square lattice (Fig. 2). Using the same methods 
as before, the expression for the GE follows. D 



5. Brief discussion of the honeycomb lattice results 

Let us comment briefly on the above results. First, we 
have seen that the GE is, in this case, different for the 
|0,0) and |+,+) states. One could also infer from here 
that the topological contribution to the GE is also differ- 
ent. However, we wish to stress that the topological GE 
is a long-range contribution, and can only be extracted 
reliably after appropriate blocking of the spins. Second, 
it is clear that the same conclusions that applied to ex- 
cited states and other ground states on the square lattice 
apply also to the honeycomb lattice, i.e. once chosen a 
ground state basis, all the states in the basis as well as the 
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FIG. 6: (color online) Mapping of the toric code from a square 
lattice S to a honeycomb lattice S (Levin- Wen string model) , 
and back to a square lattice E'. Each red arrow represents a 
CNOT operation that disentangles qubits from the system, as 
explained in Fig. 2. The arrows go from controlling to target 
qubits. 



corresponding (quasiparticle) excitations have the same 
amount of entanglement as the "reference" ground state 
(which in this case is either |0,0) or |-|-,-|-)), as long as 
they are related by local operations. 

Moreover, the blocking of spins in the honeycomb lat- 
tice produces similar results for the GE as for the square 
lattice. This is so, because in just one RG step we can 
map the honeycomb to the square lattice. Let us explain 
this in more detail. It is well known that the toric code 
in the honeycomb and square lattices can be mapped to 
one another by simple RG (disentangling) operations act- 
ing locally on the lattice. These transformations (from 
the square to the honeycomb lattice, and back to the 
square) are represented in the diagram of Fig. 6. In fact, 
these moves can also be understood in terms of Entan- 
glement Renormalization steps [ ] , where some of the 
initial qubits become disentangled from the rest after ev- 
ery step. 

Importantly, after mapping the honeycomb to the 
square lattice, the GE follows exactly the laws explained 
in the previous section, and hence it is evident that 
E^ — 1 also holds in this case, after a suitable block- 



ing of the spins. 



IV. BEYOND THE TORIC CODE 

After our detailed study on the toric code model, we 
now turn to study the GE in other topological models 
that also correspond to fixed points of RG. From now 
on, some parts of our study will not be as detailed as 
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for the toric code, and will focus mainly on some of the 
fundamental properties. Yet, in order to make everything 
as clear as possible, we need to see how some of the main 
properties of the toric code apply also to other models. 
In particular, there are three main points to consider: 



will be done in Sec. VII, and quantum double models will 
be considered subsequently in Sec. VIII. 



DOUBLE SEMION MODEL 



1. Excitations and other ground states.- As for the 
toric code, excited states and other ground states that 
can be obtained by local operations acting on some ref- 
erence ground state will share the same entanglement 
properties (c.f. Lemma 1). Therefore, in many cases we 
shall only discuss the entanglement of some "reference" 
ground state. 

2. Disentangling general stabilizer states.- Very im- 
portantly, the disentangling property is not exclusive of 
the toric code only, but it holds for any stabilizer state as 
well. More concretely, any state that is stabilized by a set 
of operators can be built from a quantum circuit acting on 
some separable initial state, where the gates of the quan- 
tum, circuit involve the basis of the stabilizer group (see 
e.g. [20]). If this basis of stabilizers can be thought of as 
acting locally on a lattice, then the quantum gates of the 
corresponding quantum circuit will also be local in the 
lattice. Moreover, there will be a set of unitary gates in 
the quantum circuit that will correspond to each one of 
the elements of the group basis, and which will take into 
account the existence of that particular basis element. 
That is: removing the action of the gates correspond- 
ing to a particular basis element effectively removes that 
stabilizer from the basis of the stabilizer group (e.g. re- 
moving a plaquette or star operator from the toric code) . 

Therefore, simply by reversing the action of such a 
quantum circuit on the stabilized state, we will always 
manage to disentangle the state by means of sufficiently 
local unitary operations. At every step, the remaining 
quantum state will be the state stabilized by the remain- 
ing set of stabilizer operators. In particular, this implies 
that there is always a MERA representation for any sta- 
bilizer state (not necessarily topological), and the toric 
code model is just a particular case. 

This property turns out to be quite important, since it 
allows us to extend straightforwardly the results for the 
GE of blocks beyond the toric code: when blocking is 
considered, simply perform disentangling operations in- 
side of each block in order to remove the maximum of 
stabilizers within each block, so that the remaining sta- 
bilizers will always be between spins in different blocks. 
By construction, the remaining state will usually obey an 
area law for the GE of blocks, plus a possible topological 
term. Because of this, from now on the GE of blocks will 
not be discussed in detail except for a couple of examples. 

3. The non-abelian case.- It is indeed possible to gen- 
eralize our results to systems with non-abelian symme- 
tries, such as quantum double models. However, in order 
to deal with these models, it will be worthwile to develop 
part of our formalism in the most general fashion. This 



The double semion model is given by the spin model 
on the honeycomb lattice 



-^DS - ^zl^'^^zl^'p' 



(24) 



where Ag and B' are mutually commuting and given by 
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(25) 



This is, the star operators are the same as for the toric 
code, whereas the plaquette operators are related via 



K 



-B 
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(26) 



with Bp the toric code plaquette operator. As for the 
toric code, star and plaquette operators satisfy the non- 
local constraint 



n^^-n^;-i 



(27) 



This model is known to be also topologically ordered, yet 
its topological phase differs from the one of the toric code 
model. More precisely, while the toric code corresponds 
to the topological phase of a Z2 gauge theory, the double 
semion model corresponds to a f/(l) x t/(l) Chern-Simons 
theory [ ]. 

One of the ground states of this double semion model 
is almost identical to the toric code ground state |-|-, -t-), 
except that each term in the decomposition Eq. (22) is 
multiplied by a factor (—1)"^^, where X^ counts the num- 
ber of loops formed by spins in the |— ) state. More specif- 
ically, this ground state is given by 



1 



;ds 



EM) 



x,^ 



(28) 



geGp 



where again Qp is the group of all possible products of 
plaquette operators B'p. Of course, the state |+,+) is 
by construction a stabilized state of the group Qp. We 
provide now the following Theorem: 

Theorem 4. For the double semion model in a hon- 
eycomb lattice E, the GE of spins for the ground state 
|+,+)j-ig is given by 



Er 



1, 



(29) 



where rip is the number of plaquettes (faces) in E. 
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Proof. Remember the ground state |+,+) for the toric 
code in the honeycomb lattice. This state can always 
be written as |+,+) = '^Cijk...\ijk...) with Cijk... > 
in the local basis {|+), |— )} for every spin. From here, 
it follows that \J2,jk..Pi<lj-Ctjk..\ < Eijfc... b»9j--Cijfc|, 
where the equality holds ii pi,qj, ... are all non-negative 
up to a common phase, and the maximum is achieved by 
choosing pi,qj, .. to be non-negative. Next, let us define 
dijk... = Cijk..e^'^, where (j> is an arbitrary phase. Then it 

follows that \Y.ijk..Pi'lj--dijk..\ < E^jfc.. \P^qJ■■■■C^jk\■ A 

maximization over the left hand side gives the maximal 
overlap, say, for the double semion model by choosing 
appropriate phase factors and coefficients. But in any 
case, this will be smaller than maximizing the right hand 
side whose maximum is, in this particular case, the max- 
imal overlap for the toric code. From here the result just 

follows. n 

From the above Theorem it follows that the topolog- 
ical contribution is E^ = 1, which is again equal to the 
topological entropy for this model. Moreover, even if we 
do not discuss it here in detail, we expect the GE of 
blocks to obey similar laws as for the toric code model, 
since the ground state can also be disentangled by lo- 
cal moves [27]. Similarly, the entanglement properties of 
other ground states and excited states obtained by acting 
locally on |-|-, -I-), remain identical to those of |+, -I-) (c.f. 
Lemma 1). 



VI. COLOR CODE MODELS 

Topological color codes [1 2] are defined as follows: 
qubits reside on vertices (not links!) of a lattice that 
is a 2-colex. A 2-colex E is a 2D lattice embedded in 
a torus of arbitrary genus g, with the following proper- 
ties: (a) every vertex of the lattice has degree 3, i.e. pre- 
cisely 3 edges meet at each vertex; (b) the faces of the 
lattice are 3-colorable, i.e. the faces can be colored with 
3-colors (e.g. red, green and blue) in such a way that no 
two adjacent faces have the same color, see e.g. Fig 7. 
Two types of commuting operators are defined on each 
plaquette (or face) p, 



B^^Yla^^^ Bl^^a 



\3\ 



j£p 



jep 



The Hamiltonian of the color code is thus 



^cc - - Z^^^P + ^P^- 



(30) 



(31) 



For this model in 2D lattices, the ground space has de- 
generacy 2''b ^ where g is the genus of the underlying Rie- 
mann surface. For instance, on the torus this degeneracy 
is 16. 

Let us define the group Gp as the group generated by 
the products of all possible plaquette operators B^ for 
all three colors. Notice that the product of all plaquette 




FIG. 7: (color online) Color code on a honeycomb (brickwall) 
lattice S, which is an example of a 2-colex. There is a spin 
on every vertex, and every plaquette is colored with either 
red, green or blue. Non-local string operators wi,W2, W3 and 
TO4 are also shown, where periodic boundary conditions on a 
torus are assumed. Grey dots indicate the sites where string 
operators are acting. 



operators B^ corresponding to the same color c corre- 
sponds to the same element of the group Gp (namely, 
the action of a ax operator on all the qubits). Therefore, 
we have the non-local constraint 
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(32) 
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The above relation implies that the size of the group is 
"p~2^ with Hp the number of plaquettes in lattice 



iXl 



The ground level subspace of the color code model 
is constructed in a very similar way as we did for the 
toric code. Specifically, for the color code on a torus the 
ground space reads C — span{\i,j,k,l), i,j,k,l — 0, 1}, 
where 



\hj,k,l) = 



1 



\G^\ 
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see,? 
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(33) 



Operators {wi^'W2tW^,wa\ above correspond to non- 
contractible loop operators, see Fig. 7.b. We do not en- 
ter here on the specific action of these operators on every 
spin in the loop, and just mention that it is local. The 16 
vectors \i,j,k,l) are orthonormal and stabilized by Gp , 
and therefore form a possible basis of the ground level 
space of the model on a torus. 

As for the toric code, the action of loop operators cor- 
responds to local unitary operations on a reference state 
|0, 0, 0, 0), and therefore do not change the entanglement 
properties (c.f. Lemma 1). Therefore, from now on we 
shall consider only the entanglement properties of the 
ground state |0, 0, 0, 0), defined as 
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In what follows we show two different approaches to 
compute the GE for this model. The first approach is 
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based, as all the calculations before, on finding upper 
and lower bounds for the GE. This approach works well 
for e.g. finding the GE of spin and blocks on regular 
lattices such as the honeycomb lattice. The second ap- 
proach uses a different technique, building on recent work 
[28]. The main idea is to exploit that color code ground 
states are instances of so-called Calderbank-Shor-Steane 
states of self-orthogonal type. This refers to the fact that 
|0,0,0, 0), when expanded in the standard basis, has the 
form 



(a) 



|0, 0,0,0) 



1 



El-) 

xes 



(35) 



where S C Zj is a collection of bit strings with a partic- 
ular structure i.e. it is a self-orthogonal classical linear 
code. This property is in fact displayed by color code 
states associated with arbitrary lattices (2-colexes) and 
for surfaces of arbitrary genus. As a result, the second 
approach will lead to a general calculation of the GE for 
arbitrary color codes. 



A. Bound method 

Let us consider the color code on a honeycomb lattice 
E. For this lattice, and with this method, we can provide 
the following Theorem: 

Theorem 5. For the color code model in a honeycomb 
lattice S, the GE of spins for the ground state |0, 0,0,0) 
is given by 



En = rir. 



2, 



(36) 



where rip is the number of plaquettes (faces) in E. 



Proof. First, since the group Qp has size \Gp\ — 2"''^^, 
we obtain a lower bound on the maximal overlap A^ax > 
l^p^l"^ ; and hence an upper bound on the GE given by 
Eg i£ Up — 2. Second, in order to derive a good lower 
bound on the GE, we consider the natural bipartition for 
the honeycomb lattice, see Fig. 8. a. One can convince 
oneself easily that any group element in Qp must have 
nontrivial support contained in both the A and B sets, 
and therefore the subgroups Gp{A) and Gp{B) are both 
trivial and contain only the identity element. For the 
color code, the reduced density matrix of A also satis- 
fies the property p\ - ilG^ {A)\\g^{B)\/\g^\)pA [29]. 
Therefore, A^^^^ < l/\g^\, and hence Eg > Up - 2. We 
therefore arrive at Eg = rip — 2, which proves the theo- 
rem, n 

Some comments are in order. First, we see that the 
term of topological origin is -E^ = 2, which is once again 
identical to the topological entanglement entropy ["_"]. 
Second, this topological geometric entanglement is twice 
that of the toric code, which is a direct consequence that 
the color code model we consider here is equivalent to 
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FIG. 8: (color online) (a) Natural bipartition of the honey- 
comb lattice E into A and B sets (respectively red and blue), 
(b) Possible blocking and bipartition of lattice S, obtained 
after applying disentanglers corresponding to the stabilizers 
within each block to the honeycomb lattice E of the color 
code. We also show the relevant bipartition into A and B 
sets (again red and blue). 



two copies of the toric code under some given RG pro- 
cedure [I-')]. Third, a similar result holds also for the 
ground state associate with the group Q^ generated by 
the product of B^ operators. Fourth, excitations and 
other ground states related by local operations will have 
the same entanglement. And fifth, the GE of blocks will 
consist again of an area-law term plus a topological con- 
tribution which will be again E-y = 2. This can be seen, 
e.g. by choosing the blocking from Fig. 8.b, and appro- 
priately disentangling qubits within each block. 



B. Color codes and self-orthogonal classical codes 

Here we provide a generalization of theorem 5 which 
will allow to compute the GE of color code states in a 
much broader setting i.e. for arbitrary lattices (2-colexes) 
and arbitrary genus. The technique used to arrive at this 
generalization is entirely different from the proof of the- 
orem 5. It involves connections between color codes and 
classical coding theory which were found in the recent 
work [::s]. 

We consider a color code defined on an arbitrary 2- 
colex embedded in a surface of arbitrary genus. In anal- 
ogy with (34) we consider the state 
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We will show here that the result in theorem 5 extends 
to all states |0): 

Theorem 6. Consider the color code ground state |0) for 
an arbitrary 2-colex E embedded in a surface of arbitrary 
genus. The GE of spins for this state is 
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(38) 
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Before proving the result, we develop some preliminary 
material. A linear subspace C C Z2 is called a (classical) 
binary linear code of length n. The elements of C are 
called its codewords. For every binary linear code C of 
length n, define an n-qubit state 



\C):= 
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We shall refer to any state of this kind as a Calderbank 
Shor-Steane (CSS) state. It easily follows from (37) that 
|0) is a CSS state. More precisely, we have |0) 
where the linear code C^ is defined as follows 



writing 



every element g & Gp as 



X"" 



(40) 



with Mi e Z2, the code Cp is simply given by the collec- 
tion of all bit strings (mi, . . . , m„) arising from elements 
of Gp in this way. In particular, one has 



\c^\^\g^\ = r 



(41) 



which will be important below. 

It turns out that color code states are CSS states of a 
particular kind. A linear code C is called self-orthogonal 
if every two codewords are orthogonal: s'^t = for every 
s,t G C. A CSS state for which the underlying classical 
code is self-orthogonal will be denoted CSS state of self- 
orthogonal type. In [ ] the following was shown: 

Lemma 3. Consider a color code associated with an ar- 



bitrary 2-colex E. Then the ground state 
CSS state of self-orthogonal type. 
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Next we provide a calculation of the GE of arbitrary 
CSS states of self-orthogonal type — by lemma 3 this will 
yield the GE of color code ground states. Note how- 
ever that not all CSS states of self-orthogonal type are 
color code ground states i.e. our results extend beyond 
the color code setting. It is also interesting to point out 
that toric code ground states, even though they are CSS 
states, are not of self-orthogonal type. 

Consider an n-qubit system and let I*!") == (^ I'/''*') be 
an n-qubit product state, where \(t>'-^') — ai\0) -\- bi\l). 
Define an associated n-qubit tensor product operator as 
follows: 



A 



I Ai , where Ai 



tti hi 



(42) 



Now consider an arbitrary self-orthogonal linear code C 
of length n and the associated n-qubit CSS state |C). 
The following lemma, proved in [ ], relates the overlap 
(CI*) to the expectation value (C\A\C): 

Lemma 4. Let C be a binary self-orthogonal linear code 
of lenghth n and let |<I>) be an n-qubit product state. Then 



(43) 



\C\-{C\^)^{C\A\C). 



Lemma 4 will now be used to prove: 

Theorem 7. Let \C) be an arbitrary n-qubit CSS state 
of self-orthogonal type. Then the CE of spins is 



EG{\C))^\og\C\. 



(44) 



Proof. As before, let Ai„ax denote the maximal value of 
|(C|<i>)| when |<I>) ranges over all product states. First, 
taking an arbitrary u £ C yields |(C|u)| = \C\~2 so that 
we obtain the lower bound Amax > |C'I~^- Second, since 
\C) has real and nonnegative amplitudes in the standard 
basis, the maximal value of |(C|$)| will be reached by a 
product state which has real, nonnegative amplitudes as 
well. Consider an arbitrary product state |<I>) with real 
and nonnegative amplitudes and the associated tensor 
product operator A as in (42). Using lemma 4 we find 



\{cm 



= \{C\A\C)\ < — ^1 

\c\ ' ' " - v4^ 
1 



^11 



nii^^ 



(45) 



where || • || denotes the operator norm. Since each \4>^'^') is 
a unit vector with real and nonnegative amplitudes, each 
Ai is a real orthogonal matrix. Hence ||yli|| ~ 1. This 
proves the upper bound Amax < |C'|^2 which matches 
the previously found lower bound. The result now follows 
immediately. D 

Combining lemma 3, theorem 7 and Eq. (41) immedi- 
ately proves theorem 6. 



VII. GENERAL FORMALISM FOR THE 
BOUND METHOD 

In order to extend our results to the non-Abelian set- 
ting (e.g. quantum double models), it is worthwhile to 
generalize the bounding technique that we have used so 
far. Looking back, we have calculated the GE of quan- 
tum states with a certain (Abelian) group symmetry, and 
our bounding strategy involved two steps. First, we chose 
suitable product states to lower bound the maximal over- 
lap (and hence to upper bound the GE). Second, we 
picked certain bipartitions and determined the maximal 
overlap with bipartite product states in order to upper 
bound the maximal overlap with multipartite product 
states (and hence to lower bound the GE). As we have 
seen, this procedure gives tight bounds in almost all the 
cases considered. 

We shall now distill the essence of this bounding strat- 
egy and derive general bounds on the GE of quantum 
states with a non-Abelian group symmetry. Other than 
placing all of our previous examples under a common 
roof these bounds may even be of independent interest 
beyond the present study of topological order. 
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A. Overview 

Let G be a finite group. Suppose its n-fold direct prod- 
uct tj = G^" acts on some m-partite Hilbert space V. = 
v.!®- ■ ■^'Hm via product operators O* = O^ (8) • • • O^ . 
Consider states of the form 



teS 



(46) 



where |<i>) — | </>'"'"') ®---® 10^™') is some reference product 
state such that any two local states 0\ |</'''^1) and O^ l^'*^') 
are either identical or orthogonal. 

We can define the global stabilizer Qi^ C ^ as the sub- 
group which leaves the reference state |<i>) invariant. Note 
that Qi^ need not be normal in general, so the collection 
Q' :— Q/Qg, of cosets is not necessarily a group. It is clear 
that the state 



^^ teg' 



(47) 



is correctly normalized. Picking any component 0*|<I>) 
we get a trivial lower bound on the maximal over- 
lap Aniax(^) and hence by (6) an upper bound EG{'i) < 
nlog2|G|-log2|^<[,|. 

Now let H — Ha ^ "Hb a bipartition of the Hilbert 
space into subsystems A and B with the reference state 
|*&) = \4'a)^\4'b) partitioned accordingly. For each oper- 
ator O* we define a truncated operator O^c as the restric- 
tion of O* to 7ix and the identity in Hx- We further- 
more define the local stabilizer Q^(X) C G for subsys- 
tem X as the subgroup whose truncated operators O^ 
leave the reference state |<i>) invariant. It is clear that 
the global stabilizer is contained in any local stabilizer, 
so G'^{X) := G<i,{X)/G<!, makes sense. Note that again 
G'^{X) is not necessarily a group. Nevertheless, for bi- 
partitions with the property 



G^iA) = G^{B) 



(48) 



we will find (see Lemma 6 below) that the largest 
eigenvalue of the reduced density operator p^ 
equals |t/$(^)| /(|t/||t/$|) which by (7) immediately 
implies a lower bound on the GE. We can summarize 
both bounds on the GE in the following 

Theorem 8 (General bound). Let |^) a state of the 
form (46). // there is a bipartition H — Ha ® Hb obey- 
ing (48), the geometric entanglement of |^) is bounded 
by 



'log. 



15*1 



<ii;G(*)-nlog2|G|<-log2|g*|. (49) 



As is evident the derivation of this general bound is 
independent of the reference state |4>) or any physical 
model in which the state |^) might arise, including any 
geometry or topology that might be involved. So the 



problem of obtaining actual bounds for specific states 
in specific models with specific geometries (topologies) 
reduces to the much simpler problem of analyzing the 
interplay between the group G, the global stabilizer fj$ 
and the local stabilizers G<i>{X) in the particular case at 
hand. Additionally, specific models will often allow us to 
characterize these stabilizers in purely geometric (com- 
binatorial) terms, which provides a constructive way to 
find bipartitions obeying (48) . Remarkably, this strategy 
even yields exact values for the geometric entanglement 
in many interesting cases as we will show (and have seen 
already). Namely, with a mild assumption in addition 
to (48) we immediately get 

Theorem 9. Let j^P) a state of the form (46). If there 
is a bipartition T-L = Ha ® Hb such that 

Gi{A)^G<s>{B) = Gi, (50) 

the geometric entanglement of |^) is given by 

£;G(*)=nl0g2|G|~l0g2|e*|. (51) 

We would like to mention that if G is Abelian then 
both G' and G'^p {X) are indeed groups and one can deter- 
mine the spectrum of pA directly in terms of these. This 
has been discussed in Ref. [6] for the special case G = Z2 . 
There our G' is denoted by G and our G'<^{A) by Ga (and 
similarly for the other subsystem B). 



B. Some Details 

Here we will supply some additional details glossed 
over before. First the precise definitions of global and 
local stabilizers. 

Definition 1 (Global stabilizer). 

G^:={teG\0'\^) = \^)}. (52) 

Definition 2 (Local stabilizer). 

G<i>{X):^{teG\0'xm^m}. (53) 



This shows how we obtain arbitrary reduced density 
operators from the state (47). 

Lemma 5 (Reduced density operator). For any biparti- 
tion H = Ha ^ Hb into subsystems A and B one has 

PA^T^.Y. E 01\<Pa){<Pa\07'' . (54) 
Proof. Avoiding the sum over cosets, we may rewrite (47) 



as 



I*) 



— =^^0^10^) ® 0*5 |0b) 
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hence the global density operator reads 
1 



P 



mQ^\ 



J2 01\cI,a){<I>a\01' '8>0-s\^b){cI>b\0^b" ■ 



u.v^Q 



By assumption {O^ \(f)B) \ t S Q /Q<i>{B)} is a set of mu- 
tually orthogonal states in subsystem B so we can use it 
to take the partial trace: 



trs(p) 



1 



1^11^*1 



^ 01\cj,a){cI^a\0T" ■ {M01\<Pb). 



u.v(EG 



Since the local expectation value is non-zero precisely for 
V G G<p{B) we obtain 






We turn to those bipartitions obeying (48). 



vu 
A 



D 



Lemma 6 (Spectrum of reduced density operators). Let 
% = Ha ® T~iB o- bipartition such that Q^{A) — Q^(B). 
Then the spectrum of pA is flat with non-zero eigenvalues 



1^11^*1 ■ 
Proof. From Lemma 5 we conclude that 



(55) 



PA 



\G\\g<i. 



J2o\\cj)A){^A\oy, 

teg 



D 



and now the observation 

is enough to prove the claim. 

VIII. QUANTUM DOUBLE MODELS 



Consider a planar graph T = {V, E, F) with vertices V, 
edges E and faces F, and sizes \V\ — rig, \E\ = n and 
|F| — Up. The Hilbert space of the quantum double 
model is defined asH '■— Hi iS) ■ ■ ■ 'S) Hiei where each 
local Hi ~ CG. Vertex projectors (acting on stars) read 



A. 



' ' sec 



where each individual vertex operator A^ acts on a ver- 
tex (star) s £ V by multiplying all edges meeting at s 
by g E G from the left, provided all edges point towards s. 
Otherwise an edge is multiplied by g^^ from the right. It 
is clear that vertex operators on different vertices com- 
mute. Similarly, plaquette projectors are defined as 



Bp = 5{ W gi,e), 



(56) 



which selects configurations where the ordered product 
of group elememts along an oriented circuit Cp around p 
is the unit element of G. The Hamiltonian is the sum of 
these mutually commuting projectors over vertices and 
plaquettes, namely 



H 



D(G) 



sev 



A.. 



E 

peF 



Br, 



(57) 



We identify the direct product group from our general 
discussion in the previous section as G :— G^'^' and ob- 
tain its natural action on the Hilbert space by collecting 
individual vertex operators into the joint vertex operator 

o* :- n ^'° 

sev 

where t ~ (gi, ■ ■ ■ , g\v\) G G- Clearly, any joint vertex 
operator has product form: given a directed edge (s, s') 
with value \x) an element t = {. . . ,gs, ■ . . ,gs', • . . ) S ^ is 
easily seen to act locally as 



^(s,s') \^) 



\gs'xg; 



(58) 



thus O* = Ol 



)0 



\E\- 



1^1 

We can obtain a ground state of the quantum dou- 
ble model by choosing the reference product state |e) :— 
{^^|e) and projecting it on the common -|-1 eigenspace 
of all vertex operators 



i")=n^^i^)=^Eo 



sev 



\Q\ 



(59) 



teg 



i along Cp 



In order to proceed we need to figure out the actual 
form of the global and local stabilizers in the quantum 
double model for our particular reference state. As for 
elements of the global stabilizer Qe, it follows immedi- 
ately from (58) that ga' — gs iff vertices s and s' are 
connected by an edge. Hence Qe only depends on the 
connectivity of the graph F. Note that Q^ is not normal 
in Q generally. As far as the local stabilizers Qe{X) are 
concerned, we can equivalently define them combinatori- 
ally. Partition the vertices V into clusters Vi{X) based 
on whether vertices are connected by an edge in Ex. 
(That is, two vertices are in the same cluster iff they can 
be connected by a path which lies completely in Ex.) 
Then t = (gi, . . . ,g\y\) G Qe{X) iff its components gs are 
constant on clusters. It turns out that local and global 
stabilizers are related quite favourably: 

Lemma 7 (Joint local stabilizer). Ge{A) D Ge{B) = Ge- 

Proof. We want to cover the vertices V with a 
set {^^(X)} of overlapping clusters obtained from any 
of the subsystems A or B. Then if t = {gi, . . . ,g\v\) S 
Ge{A) n Ge{B) its components gs must be constant on 
each cluster, and hence constant on the whole of V be- 
cause of the overlaps. Clearly, we can obtain a global 
cover by constructing a local cover for the vertices along 
each (edge) path 7. 
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So let 7 be a path starting at vertex s and let Vi (A) be 
the enveloping cluster of s as given by subsystem A. If 7 
never leaves this cluster we are done. Otherwise we may 
assume that s lies at the boundary of Vi{A), hence we 
will reach a distinct cluster V2{A) by advancing a single 
edge (s, s') along 7. Clearly, this edge must be in Eb and 
thus there exists a cluster Vi{B) containing both s and s' . 
Then by induction {Vi{A),Vi{B), V2iA), . . . ,Vi{A)} is a 
local overlapping cover along 7. D 

So in order to find a suitable bipartition with GeiA) — 
Ge{B) = Qe we simply need to select a subset Ea of 
edges which connects all vertices into a single cluster. It 
is then clear that \Qf.\ = \G\. Indeed, for the square and 
triangular lattice such bipartitions exist, see Fig. 3(a) and 
Fig. 5 for examples. Using Theorem 9 we then obtain 

Theorem 10. For any quantum double model on a 
square or triangular lattice E, the ground state jfi) has 
the geometric entanglement 



EG = a\V\-\og^\G\. 



(60) 



where a = log2|G|. 



Note that \V\ = rig (number of stars) is proportional 
to the volume of the surface. 

From [ ] and the proof of Theorem 2 it becomes 
clear that we can calculate the GE of blocks of linear 
size fc > 2 directly from the renormalized graphs Tk 
shown in Fig. 2.e. On a square lattice of linear size 2k X 
there are N^ = 2A^ such blocks, and the renormalized 
graph Ffc has \Vk\ = Nk{k + 1) vertices. We merely state 
the result for the GE of blocks. 

Theorem 11. For any quantum double model on a 
square lattice S, the ground state \fl) has the geometric 
entanglement of blocks of size k 



EG = a\Vk\-log2\G\ 



(61) 



whe 



log2|G|. 



Clearly, \Vk\ is proportional to both the renormalized 
volume of the surface and the boundary area of a single 
block (with a natural value of 4A:). This is exactly the 
generalization of Eq.(14) to the non-Abelian case. As a 
final remark, we notice that the topological contribution 
to the geometric entanglement for quantum double mod- 
els is given by E^ — log2 \G\. This again coincides with 
the corresponding value of the topological entanglement 
entropy. 



a perturbation? In this section we address briefiy this 
question, and arrive to the conclusion that, indeed, E^ is 
a robust property of the topological phase. 

There are several ways of checking the robustness of E^ 
under perturbations. One possibility is to perform a nu- 
merical analysis. We shall not pursue this approach here, 
even though the procedure for this numerical evaluation 
is already clear (e.g. by finite-size scaling of the GE of 
blocks in cylindric geometries) and will be the subject of 
a future work. Another option is a perturbation theory 
analysis. Yet, a more intuitive alternative is the following 
argumentation based on RG fixed points: all the models 
considered here are RG fixed points and hence represen- 
tatives of their respective topological phases. The long- 
distance properties of any state in one of these phases do 
not change under local RG transformations and, thus, are 
equivalent to those of the fixed point. This, in particu- 
lar, is true for the long-range pattern of entanglement, 
and hence for E^. Thus, any non-relevant and short- 
range perturbation driving the Hamiltonian away from 
the fixed point will produce ground states with the same 
E^. Nevertheless, we expect a change in the short-range 
pattern of entanglement , and hence in the non-universal 
term corresponding to a boundary law. For short-range 
perturbations the change involves modifications of the 
prefactor of the boundary law as well as the possible 
appearance of subleading 0{L~'^ ) corrections. Hence, 
Eq.(4) applies when away from the fixed point. 

To double-check the above claim, we now perform a 
simple perturbation theory analysis of the robustness of 
the GE of the toric code model under external short- 
range perturbations such as magnetic fields. This anal- 
ysis provides upper and lower bounds for the GE, and 
complements the above argumentation on the robustness 
of E^ based on RG. 

Let us then add a perturbation to the toric code Hamil- 
tonian on the square lattice, and see how the ground state 
changes. For simplicity, we consider the case of an infinite 
plane, where the ground state |0,0) is non-degenerate, 
and hence we can use non-degenerate perturbation the- 
ory. The perturbed Hamiltonian will be 



H^ = H' 



TC 



AK 



(62) 



where i?TC is the toric code Hamiltonian, V is the pertur- 
bation, and A ^ 1. Non-degenerate perturbation theory 
says that the new ground state can be approximated as 



IX. E^ AWAY FROM FIXED RG POINTS 

So far we just discussed the topological contribution 
to the GE for specific ground states of models that turn 
out to be RG fixed points, and therefore representative 
of their respective topological phases. But, what if we 
are away from the RG fixed point? Is E^ robust under 



|0,0)^»|0,0) + Ai:%MlM|,,,,„,0) ,63, 



£-0,0 — E, 



(p.c 



where Eqa is the ground state energy and E^j^c is the 
energy of the excited state |0, c, 0, 0) . 

We now consider the case in which the perturbation is 
an homogeneous magnetic field, e.g. in the x direction. 



V - 






(64) 



(the case of z and y directions can be considered simi- 
larly) . It is easy to check that in this case, the normalized 
perturbed ground state becomes 



|0,0)'«C |0,0)--^aW|0,0) 



(65) 



j=i 
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with A the energy gap to create a pair of flux and anti- 
flux quasiparticles, and C = (l-|-riA^/A^)~^/^ a normal- 
ization constant. 



Our aim now is to estimate the maximum overlap of the 
previous state with a product state of blocks of boundary 
size L. This can be done as follows: first, and as in 
the unperturbed case, we apply CNOT operations locally 
inside of the blocks so that qubits are disentangled in the 
unperturbed ground state. By doing this, we can focus 
on the entanglement of the state 



|0,0)d'i 



disentangled 



C |0,0)--^a[^-l|0,0)U|ei) 



(66) 



j=i 



where jcfc) is the quantum state for the fc-th disentangled qubit. The above equation is indeed equivalent to 



|0,0) 



A 
discntanglci 



d«C |O,O)0|ei) 



£5[-'l|0,0)®|ei) 



j=i 



|ep) - ^|0,0)|a;i,...^p/ , , 



(67) 



where Sx is the total spin in the x direction for the L the fact that both bounds leave the topological compo- 
spins in the boundary of block j , and 



nent E^ 



1 untouched seems to indicate that this is 



wi, 



E-. 



m 



ei 



® |ep) 



(68) 



Now we find upper and lower bounds to the maximum 
overlap of this state with a product state of the blocks. A 
lower bound can be easily obtained by the product state 
|0^®("-p) (g, |e^) (g) . . . (g) \e.p). Noticing that |efc) is either 
|0) or |-|-) [24], and that the |-t-) contributions come only 
from qubits close to the boundary of the block, we have 
that 



qe^r^/^ (l - '^) < A£--I , (69) 

for some positive w = 0(1) constant. The following up- 
per bound can also be found easily: 



A^i-'^^i < c\g\-' 



-1/2 



1 



rn,L\ 

A 



(70) 



Using the above bounds, one can check that for the GE 
we obtain, in the limit A ^ A and L':^ \, 



2|^)n,L-l>i^^>Q-|)n,L- 1.(71) 



The above equation is compatible with a leading change 
in the GE in the prefactor of the boundary law. Also, 



actually robust under the perturbation. Moreover, im- 
plementing finite-L corrections to these bounds provides 
0{L^'^ ) corrections, which is consistent with our previ- 
ous claims. 



X. CONCLUSIONS 

In this paper, we have studied the GE in various 
topologically-ordered states that correspond to fixed- 
points of RG. These are the toric code, double semion, 
color code, and quantum double models, in a variety of 
lattices. Generically, we found that the GE is typically 
composed of a boundary law term times the number of 
blocks plus a topological term. This is remarkable, since 
it is a signature of topological order in multipartite en- 
tanglement properties of the state, rather than bipartite. 
The topological term for the color code is twice that for 
the toric code, which is consistent with the recent re- 
sult that the former model is equivalent to two copies of 
the latter [l.'i]. Away from RG fixed-points, we argued 
that the same type of behavior holds up to a possible 
subleading term, and thus the topological contribution is 
a robust property of topological phases of matter. The 
numerical evaluation of E^ for non-analytic cases is pos- 
sible. This could be done e.g. in the context of Tensor 
Network methods (such as PEPS algorithms) in cylindric 
geometries of increasing cylinder size, as is usually done 
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for the topological entropy [■*"]. We leave this also for 
future work. 

In all the cases considered here, the topological GE 
turns out to be equivalent to the topological entropy. It 
would be good to prove, in general, whether this is al- 
ways true or not, in order to understand if the GE can 
provide more information about topological order than 
the one that is already in the topological entropy. How- 
ever, even if this were not the case, the GE would still be 
a useful tool to extract the quantum dimensions of the 
anyon model associated to a given topologically ordered 
state. For scenarios where reduced density matrices are 
hard to evaluate numerically, one may thus believe that 
looking at the GE may be more efficient. 

Finally, we believe that other models could possibly 
be analyzed with the techniques that we used in this pa- 
per. Specifically, similar results should also apply to the 
A-phase of Kitaev's honeycomb model [ ] for which the 
toric code is (in some limiting cases) an effective model, 
as well as spin-liquid states with an emergent Z2 gauge 
symmetry [.32]. The possibility of studying the ground 
states of string-net models of Levin and Wen [ ] in gen- 
eral, as well as topological quantum field theories [ ], is 
also left for future investigation. 
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Appendix A: Relative entropy of entanglement 

We shall see in the following that for most of the 
ground states discussed in this paper, they can also be 
characterized by the so-called relative entropy of entan- 
glement and in fact their value is identical to that of GE. 
The relative entropy S'(p||cr) between two states p and a 
is defined via 



which is evidently not symmetric under exchange of p 
and a, and is non-negative, i.e., S{p\\a) > 0. The rela- 
tive entropy of entanglement (RE) for a mixed state p is 
defined to be the minimal relative entropy of p over the 
set of separable mixed states [33, 34]: 

Er{p) EEminS'(p||cr) = min Tr (p log2 p - p logj ct) , 

(A2) 
where V denotes the set of all separable states. With 
any separable state a, we can obtain an immediate upper 
bound on Efi{p) < Tr(plog2 p— plogj cr). For pure state 
p = |vl/)(vl/|, such as the ground state of the toric code 
and the color code, we have further that 



i^fl(|*)(*|)<-(*|log2a|v|/) 



(A3) 



In general, the task of finding the RE for arbitrary 
states p involves a minimization over all separable states, 
and this renders the computation of the RE very diffi- 
cult. It has been shown that for pure state |^), the its 
GE lower bounds its REE, i.e., ErI^) > Eci'i') [35]. 
Using the technique in Ref. [ "], we show that the GE 
is identical to the REE for the ground states considered 
here. 

Let us illustrate this by the |0,0) ground state of 
the toric code on the square lattice. Since |0, 0) = 
J2geg, ff|0...0)/-\/|t/s|, where Gs is the group generated by 
star operators, each component \g) = g|0...0) is a product 
state. We can thus construct a mixed separable state 



1 

Ws 






(A4) 



and hence obtain an upper bound on the REE, 



S{p\\cr) = Tr(plog2p-plog2cr) 



(Al) 



i?fl(|0,0)(0,0|) < log2|t;,|(0,0|E (l.9>(5l)|0,0) 

gee. 
= log^\Gs\- (A5) 

This upper bound turns out to be identical to the value 
of the GE, a lower bound on REE. Therefore, we have 
that Sfl(|0,0)(0,0|) = ^G(|0,0)(0,0|) = n,-l. Such an 
equality of GE and REE can also be understood from the 
group symmetry [30]. We can straightforwardly apply 
similar arguments to ground states of the double semion, 
color code and the quantum double models. 
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